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SUMMARY 


This paper presents a discussion of the concept of saodal densit\ 
or density of eigenvalues of various structural elements of engineering 
importance. Expressions and graphs are presented that can be used to 
estimate the average modal densities of these elements and are valid 
for elements having any prescribed boundary conditions. The expres- 
sions for modal density and their graphical representation were pre- 
pared from the information available in the literature, but supplemen- 
tary data were generated where required. 

Cases are considered for rods, beams, solid rectangular and 
circular plates, thin cylindrical, spherical and conical shells, composite 
structures, shallow sandwich shells, orthrotropic plates, pretwisted 
plates, plates subject to in-plane forces and shells on an elastic 
foundation. For each of the elements up to a composite structure, graphs 
are plotted using dimensionless parameters to generalize the applications 
of the results; however, for the rest of the elements, graphs are plotted 
by choosing some arbitrary dimensions to illustrate the effect on modal 
density. 
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1. INTRODUCTION 


Any continuous structure possesses an infinite number of natural 
modes of vibration and to obtain the information concerning the response 
of such a structure^ it is necessary to express the normal modes in 
a series form. However during the past few years^ there has been an 
effort to develop a new approach to these multimodal vibration problems 
that avoids the problem of expanding the response in terms of the mode 
shapes. In this approach^ sometimes referred to as "statistical energy 
analysis"^ average response levels in various frequency intervals are 
estimated without the apparent knowledge of the mode shapes and 
resonance frequencies. Instead^ what is required is a knowledge of the 
type and number of structural vibration modes occurring in a given 
frequency interval. This quantity^ the number of modes per unit 
frequency^ is called the 'modal density' of the structure. Thus the 
modal density of a structure is essentially the density of the modes of 
vibration with respect to frequency. It is an indication of the spac- 
ing of the natural modes in the frequency domain. 

When dealing with the structures excited in a very complex^ or 
random fashion it is often not only useful but necessary to resort to 
statistical energy analysis to determine the response of the structure 
to such loading. In order to apply this type of analysis it is found 
that the modal density of the structure in question must be known. 
Moreover statistical energy analysis shows promise of becoming a useful 
tool for estimating average response levels of multimodal structural 
vibrations as the modal density of a structure is relatively independent 



of the boundary conditions^ Hence in order to apply a statistical type 
of analysis to a structural response problem^ it is necessary to know 
the modal density of the basic structural elements such as rods^ beams^ 
and shells. It is therefore the purpose of this paper to discuss in a 
systematic manner the problem of the modal density in vibration prob- 
lems of some basic structural elements like rods^ beams^ plates and 
thin cylindrical^ spherical and conical shells, composite structures 
and certain shallow structural elements and present the expressions and 
graphs that can be used to estimate the average modal densities of 
these elements. 

The determination of the modal density is essentially a mathe- 
matical problem. It involves the determination of the frequency equa- 
tion for the structure under consideration from the appropriate 
equation of motion and then the summation of the resonant frequencies 
over all possible modes of vibration. This yields an expression for 
the number of resonant modes in terms of frequency. Differentiation of 
this expression with respect to frequency will then yield the expres- 
sion for the modal density in terms of the frequency. The k-space 
integration technique introduced by Courant and Hilbert (1953) is 
utilized to evaluate the number of resonant frequencies. 

In Chapter 3, the modal density for longitudinal and torsional 
vibrations of circular rods having uniform cross section is discussed 
and results are compared. 

Chapter 4 deals with the beams having constant geometry and 
properties. The expression for modal density for transverse vibrations 
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of beams is presented. The graph of normalized modal density versus 
dimensionless frequency is plotted. 

For Chapter 5y modal density expressions are presented in flat 
rectangular and circular plates and the results are discussed. 

In Chapter 6^ thin cylindrical shell is considered. The modal 
density expressions are developed for cylindrical shells following 
three different approaches. These integral expressions are then 
evaluated and plotted in dimensionless form. 

Chapter 7 deals with spherical shells. The modal density expres- 
sion is developed as a function of dimensionless frequency and a graph 
is plotted to illustrate the variation of modal density above and below 
the ring frequency. 

In Chapter 8^ expressions for modal density of thin conical shells 
are obtained based on two separate frequency equations and are normal- 
ized with respect to cone geometry and presented for the frequency 
range below the lower ring frequency and above the upper ring 
frequency of the cone. 

In Chapter 9^ the additive property of modal density for composite 
structures is verified analytically by considering an L-shaped frame 
consisting of two beams joined at right angles and the graph is plotted 
to illustrate the variation of modal density of the composite structure 
with respect to frequency.^ 

^Hart, F. D. and V. D. Desai. 1967. Additive properties of modal 
density for composite structures. Presented at the 74th Meeting of the 
Acoustical Society of America, Miami, Florida, Paper No. DD 11. 
Department of Mechanical and Aerospace Engineering, North Carolina 
State University at Raleigh, N. C. 
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In Chapter 10, modal density expressions are presented for shallow 
sandwich shells_, orthrotropic plates^ pretwisted plates^ plates subject 
to in-plane forces^ and shells on an elastic foundation. Graphs are 
plotted to illustrate the effect on modal density. Results obtained 

are discussed in detail and compared with some of the basic elements. 
Chapter 11 presents a summary of results and conclusions. 
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2. REVIEW OF LITERATURE 


The problem of determining the modal density of any given struc- 
ture is equivalent to ascertaining the distribution of eigenvalues of 
large order corresponding to high mode numbers □ A general discussion 
of the asymptotic distribution of eigenvalues for various classes of 
differential equations is given by Courant and Hilbert (1953) . 
Expressions for the number of eigenvalues up to a given bound are given 
for differential equations with one, two and three independent space 
variables. Although the treatment of the subject by Courant and 
Hilbert is approached from a basic mathematical point of view^ the 
results have direct physical interpretation. It is indicated that 
boundary conditions have no effect on the asymptotic distribution of 
the eigenvalues, 

Bolotin (1962) has also given considerable attention to the 
asymptotic method in his studies of eigenvalue determination. In 
1962 Bolotin presented a discussion of the asymptotic behavior of the 
eigenvalues for a generalized rectangular region of arbitrary dimen- 
sions, He applied this technique to the problem of plates and shells, 
where the number of eigenvalues correspond to the number of natural 
frequencies of vibration. Correction factors were also introduced to 
extend the work of Courant and Hilbert (1953) to low mode numbers where 
boundary conditions must sometimes be considered, Bolotin (1960) 
presented a detailed discussion of the effect of edge conditions on the 
vibrational modes of elastic shells. 
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In 1963 Bolotin presented a general treatment of the eigenvalue 
density problem for a general thin elastic shell of revolution with 
constant thickness in orthogonal curvilinear co-ordinates coinciding 
with the curvature lines» Bolotin again used the asymptotic method 
discussed by Courant and Hilbert (1953) in his work and obtained 
expressions for the number of natural frequencies and the modal density 
of a general elastic shell of revolution through elliptic integrals- 
The results of this work were also extended to the specific cases of 
the spherical shell and the circular cylindrical shell. Bolotin (1965) 
presented a discussion which was essentially an extension of his 
previous work in which he discussed the concentration points of natural 
modes^ as well as the effects of shear and rotary inertia. 

Without apparent knowledge of Bolotinas work^ Heckl (1962) 
developed an expression for the natural frequencies of a cylindrical 
shell using impedance methods. He then represented the number of 
natural modes by a finite sura over all possible modes of vibration 
possible up to some upper frequency. He then replaced the summation by 
an integral and obtained an approximate expression for the modal 
density of thin cylindrical shells. Heckl also presented some experi- 
mental findings in his report. 

In 1965 Smith and Lyon introduced the concept of modal density and 
discussed its application with regard to structural vibration. The 
cases of simply supported beams^ clamped bearas^ simply supported 
rectangular plates and clamped circular plates were considered in 
particular. 


6 


Ungar (1966) discussed the concept of modal density and its 
application to composite systems. He also presented a list of expres- 
sions for the modal density of some simple elastic systems of engineer- 
ing value. 

In 1967 Hart and Desai presented a discussion of the additive 
property of modal density for composite structures and verified 
analytically the validity of the additive property by considering the 
composite structure consisting of two beams joined at right angles to 
form an L-shaped frame. 

Miller and Hart (1967) made a combined analytical and experimental 
study about the modal densities of a thin cylindrical shell. Expres- 
sions for modal density were presented in integral form using three 
different methods and the validity of the results was discussed in 
detail. 

In 1968 Wilkinson presented the expressions for the modal densi- 
ties for transverse vibrations of two dimensional structural elements 
which included shallow sandwich shells, orthrotropic plates, shells on 
elastic foundation, pretwisted plates and plates subject to in-plane 
forces. The effect on the modal densities of these elements was 
illustrated graphically. 

In 1969 Miller presented the expressions for modal densities of 
conical shells based on the two separate frequency equations and 
applicable to a wide range of cone geometries and valid over a fre- 
quency range sufficiently wide to be of engineering value. Miller 
also presented the findings for the modal densities of conical shells 
obtained by experimental study. 
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Erickson (1969) presents expressions to estimate the average modal 
densities of sandwich beams and flat or cylindrically curved sandwich 
panels^ The effect of transverse shear flexibility^ orthrotropic 
shear moduli of the core, face bending stiffness^ rotary inertia and 
panel curvature on modal density is illustrated graphically over the 
wide-frequency range. Modal densities of flat rectangular sandwich 
panels having orthrotropic cores are determined experimentally. 
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3. MODAL DENSITY OF RODS 


3. 1 Introduction 

In this chapter a circular rod having the uniform cross section 
with both the ends fixed 'is considered. A rod can execute longi- 
tudinal^ torsional or transverse vibrations either individually or in 
combination. The expressions for modal density are derived considering 
longitudinal and torsional vibrations individually. 


3.2 Longitudinal Vibration 

The governing equation of motion for longitudinal vibrations of a 
rod is given by 


d^w(x^ t) _ ^^c t) 


bt 


bx 


(3,1) 


where 

E = Young's modulus 
p = density of material 
g^ = gravitational constant 

w = longitudinal displacement of a section. 

This equation is based on the following assumptions: 

1. The rod has a uniform cross section. 

2. During the vibratory motion^ the cross section normal to the axes 
of the bar remain plane and normal to the axis. 

3. The particle in a noruial cross section moves in the axial direction 
of the bar. 
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Letting = C ^ equation (3. 1) reduces to 

P L 


S^w(x,t) _ „2 5^w(x, t) 






(3.2) 


where C is the longitudinal velocity of wave propagation along the 

JhJ 

length of the rod. 

Assuming that the solution of equation (3.2) is 


w(x,t) = X(x) sinoDt y 


(3.3) 


and substituting (3.3) into (3.2) and simplifying^ gives 


iiiW * . 0 

dx 


(3.4) 


where 




(jj = frequency of vibration 

X(x) = the shape of normal mode of vibration. 


The solution fo (3.4) is given as 


X(x) = A cos \x + B sin Xx . 


(3.5) 


For a rod fixed at both the ends^ the boundary conditions are 



w(0, t) = w(-t, t) = 0 

(3.6) 

or 

X(0) = X(^) = 0 . 

(3.7) 
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Imposing the boundary conditions in (3.7)^ the frequency equation 


for this case can be written as 


miTC 

(U = ^ — m = 2^ 3^ 

where ^ is the length of the rod. 


(3.8) 


I 


(a) 


H 

1 

H ? 

— 


(b) 


Figure 3.1 Rod with fixed ends and k-space 


The wave number k^ may be defined as 


u _ nm 

- T 


(3.9) 


Hence the change in the wave number from one mode to the next is 
given by 


Ak. = 


1 I 


(3.10) 
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Since the waves in the case of a rod are propagated only along the 
length of the rod^ the k-space is one dimensional and the equation for 
the number of resonant frequency becomes 


N(od) 



(3.11) 


This gives 

N(cu) = ^ , 


(3.12) 


but equations (3.8) and (3.9) give 



(3.13) 


Therefore the expression for the number of resonant frequencies is 
N(U)) = J ( 7 ^) . (3.14) 

Defining a dimensionless frequency v as 


equation 



(3. 14) reduces 

»(v) . i . 


to 


Differentiating (3.16) with respect to v gives 


(3.15) 


(3.16) 


n(v) 


TT 


(3.17) 


This is an expression for modal density for the longitudinal 
vibration of a rod fixed at both ends and it can be shown to be appli- 
cable for arbitrary end conditions. 
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3.3 Torsional Vibrations 


The governing differential equation of motion for torsional 
vibrations of a rod is given as 


(x^ t) _ 


3 *(x,t) 


(3.18) 


where 

G = modulus of rigidity 

p = density of material 

g = gravitational constant 
c 

^ = angular displacement of the section. 


This equation is also based on the assumption of equation 
Defining 



(3.1). 


(3.19) 


where is the torsional wave velocity equation (3. 18) becomes 


s * (x, t) ^ ^2 a » (x, t) 


at 


dx 


(3.20) 


Assuming the solution of equation (3. 19) is 


0 (x^ t) = 0(x) sincut 


(3.21) 


substitution of (3.21) into (3.20) and further simplification gives 


t x^ecx) = 0 


(3.22) 


dx 
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where 



U) = frequency of vibration 

0(x) = the shape of normal modes of vibration. 

Again for the fixed rod^ the boundary conditions are 

0(0) = G(^) = 0 . (3.23) 


Applying the boundary conditions to (3.22) gives the following 
frequency equation: 

mnC 

(JO = — m - 1, 2, 3 ... . (3.24) 


Hence as in the previous section^ the number of resonant 
frequencies obtained by k-space integration is given as 


N(uj) 


= - (— ) 
TT 


(3.25) 


where 


C 


2 

T 




- ^ ] 
2(l+u) ^ 


(3.26) 


or 


■'T 


[2(l+n>] 


Substitution of equation (3.27) into (3.25) gives 
N(U)) = 

trc^ 


(3.27) 


(3.28) 
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Now introducing the dimensionless frequency^ equation (3.28) 


reduces to 

»<V) . 

TT 

Differentiating (3.29) with respect to v gives 

„(v) . . 

TT 

This gives an expression for modal density for the torsional 
vibrations of a rod fixed at both the ends. 


3.4 Discussion 

The expressions developed for modal density of rods for both 
longitudinal and torsional vibrations show that modal density of a rod 
is constant and is independent of geometry of cross sections in 
dimensionless form. 

Moreover^ the modal density of the rod for torsional vibrations is 
about 1.5 times that for longitudinal vibrations. 

Expressions generated by considering different boundary conditions 
give the same answer and hence it is independent of boundary conditions 
also. 
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4« MODAL DENSITY OF BEAMS 


4. 1 Introduction 

In this chapter^ the modal density for transverse vibration of 
beams having constant geometry and constant properties are discussed. 
The simply supported beam is considered for deriving the expression. 

The problem of simply supported beams was discussed in the 
literature. It is reproduced here to illustrate an exact way for 
developing the expression. Modal density is expressed as a function of 
dimensionless frequency and the graph is plotted in terms of dimension- 
less parameters. 


4,2 Simply Supported Beams 

The governing equation of motion for transverse vibrations of a 
beam is given by the differential equation 

El + EA ^ Q 

Sc at^ 

where 

E = Young's modulus of elasticity 

I = moment of inertia of cross section 

p = density of the material 

g^ = gravitational constant 

A = area of the cross section 

V(x^t) = deflection of the beam at any section. 

The equation (4. 1) is based on the following assumptions: 

1, The effect of rotary inertia is neglected. 

2. Shear displacement due to a vibratory force is negligible. 
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3, Cross sections are plane before strain and remain plane after 
strain^ 

4. Beam is slender. 

Assuming that the solution of equation (4.1) is 

V(x^t) - X(x) sinoot ^ (4.2) 


substitution of equation (4.2) into (4.1) and further simplification 
gives 

l^-x"^X(x)=0 (4.3) 

dx 


where 


/ 2 , 

4 _ ptu A 

u) = frequency of vibration 

X(x) = the shape of the normal modes of vibration. 


The general solution fo the differential equation (4o3) is given 
as 


X(x) = A sin Xx + B cos Xx + C sinh Xx + D cosh Xx (4.4) 

where A^ C and D are arbitrary constants. 

For the beam under consideration, the boundary conditions are 


v(0,t) 

= v"(0,t) = 0 

(4.5) 

V(t,t) 

= t) = 0 . 

(4.6) 
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The above conditions can be written as 


X(0) 

11 

>< 

O 

II 

O 

(4.7) 

x(-t) 

II 

>< 

I! 

O 

(4.8) 


Applying the boundary conditions (4.7) and (4.8) to equation (4.4)^ 
the frequency equation for the beam can be expressed as (Smith and 
Lyon^ 1965)^ 

U) = KC^ (‘^•9) 


where 


I 

K 


= the longitudinal velocity of wave propagation in the 
beam material along the length of beam 
= length of beam 

= radius of gyration of the cross section. 


Let the wave number be defined as 

(4.10) 

Therefore the change in the wave number from one mode to the 
next is given as 


Ak^ = J (4.11) 

Again in case of a beam^ waves are propagated only along the 
lengthy hence the k-space is one dimensional and the equation for the 
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number of resonant frequencies is 


N(od) 




(4.12) 


This gives 


N((JD) 



Combining equations (4.9) and (4.10) gives 



(4.13) 


(4. 14) 


Therefore the number of resonant frequencies is given as 


N(0)) 


I 

TT 



(4. 15) 



H ^ h / 

I 1 1 1 h— 

(b) 


Figure 4. 1 Simply supported beam and k-space 
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Let the dimensionless frequency V be defined as 


V = 




(4.16) 


Equation (4. 15) reduces to 


N(V) = i 


(4.17) 


Differentiating the above expression with respect to v gives 


n(v) = L . 

^ 2n ^Kv 


(4.18) 


This is the expression for the modal density for the transverse 
vibration of beams in terms of dimensionless frequency. 


4.3 Graphical Results and Discussion 
The result of the analytical development in the preceding section 
is represented graphically in Figure (4.2). From the graph it is seen 
that modal density for a beam decreases as the dimensionless frequency 
V increases and asymptotically approaches zero value as v becomes 
large. 
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5. MODAL DENSITY OF PLATES 


5, 1 Introduction 

In this chapter^ modal density expressions are developed for the 
free vibrations of rectangular and circular plates. To simplify the 
derivation^ a rectangular plate with simply supported edges and a 
circular plate with clamped edge are considered. The expressions 
obtained also hold good for other boundary conditions (Bolotin^ 1960) , 
Both the cases discussed in this chapter were readily available in 
the literature and the information was gathered for systematic repre- 
sentation. The differential equations governing free vibrations of 
plates are obtained by modifying the equations describing the static 
equilibrium to account for the inertia forces introduced. 


5,2 Rectangular Plates 

The governing equation of motion for free vibration of a 
rectangular plate is given as 


D[l-^ + 2 


o w o w- | ph 




= 0 


(5.1) 


where 

Eh3 

D = ^ = the flexural rigidity 

) 

p = density of the material 

= gravitational constant 
wCx^y^t) = displacement normal to the x-y plane 
E = Young’s modulus 
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= Poisson *s ratio 


M* 

h = the plate thickness. 

Equation (5.1) is based on the following assumptions: 

1. Rotary inertia is neglected. 

2. Cross sections are plane before strain and remain plane after 
strain. 

3. The thickness of the plate is small as compared to the other 
dimensions . 

4. No strain is suffered by the middle surface. 

5. Deflections are small in relation to the plate thickness. 

In general^ when a plate vibrates_, there are an infinite number of 
natural f requencies^, and each of them has a specific mode or shape of 
vibration associated with it. These modes are called normal modes or 
principal modes. 

Let the solution of equation (5.1) be assumed as 

w(x,y,t) = W(x,y)e^'“^ . (5.2) 

Substituting equation (5.2) into (5.1) and simplifying gives 

D7^W(x,y) - k^W(x,y) =0 (5.3) 


where 


phu) 


Considering the rectangular plate of dimensions and as shown 
in Figure 5.1, the boundary conditions for the plate under consideration 
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Let 


W(x,y) = sin sin 


(5.8) 


where m and n are integers. 

This function satisfies the boundary conditions for a simply 
supported rectangular plate. 

Substitution of equation (5.8) into (5.3) gives 


,2 .mTT.4 .nTT. 2 .nn.4 

‘'mn = ^ 2(;j-) (^) + 


(5.9) 


Equation (5.9) can also be expressed as 


mn "^2 


( 5 . 10 ) 


Hence the frequency equation for the simply supported rectangular 
plate is given as [Smith and Lyon, 1965] 


2 2 2 2 

U, = (5L" + !L^) KC 

mn .z .z L 


m = n = 1, 1, 3, . . , 


(5.11) 


where 


and are the length and width of the plate 


K = (\f 


12(l-n^) 


-) h « 0,289h = radius of gyration fpr the 
plate 


./I 


C = y — ^ = longitudinal wave velocity . 

L p 
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Defining the wave numbers and k^ as 
1 miT j T 

"i -z; “">* 


(5.12) 


the changes in the two wave numbers from one mode to the next is given 
as 


Aki = ^ and 



(5.13) 


Since the waves in the case of a rectangular plate are propagated 
along the length and width of the plate^ the k-space is two dimensionalo 
The expression for the number of resonant frequencies can be 
expressed as 

”<”> ■ Aiq^- /i" 

s 


Cylindrical coordinates can be utilized to integrate over the 
surface of the k.-space = 

Letting 

k^ = r cos 0 
k^ = r sin 0 

equation (5.14) takes the form 
I I r tt/ 2 

N(U) ) = r r rdedr 

TT 0 0 

Writing the equation (5.15) as 

1. ;p/ 2 I 

N(tu) = — — J [ J rdr ] 

0 0 


de 


(5.15) 


(5.16) 


26 



and carrying out the integration with respect to r first^ and then with 
respect to 0^ gives 


^it 2 

N(«,) = ^ r 


(5.17) 


Combining (5.11) and (5.12) gives 


U) 

KC. 


(5.18) 


Substitution of equation (5.18) into (5.17) gives 


N(u)) = 


I I 

12 0 ) 


4tt Ka 


(5.19) 


Define v = — — , then equation (5.19) reduces to 


N(V) - 


4ttK 


(5.20) 


where v is a dimensionless frequency parameter. 

Differentiating (5.20) with respect to v gives 


n(v) 


^ dN(v) _ ^^2 

dv 4^K 


(5.21) 


This is an expression for the modal density of a rectangular 
plate. 

If the plate thickness is h_, the radius of gyration is h/\[TT 
and the expression is given as 

^2V3 

“ ZtST. (5.22) 
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5»3 Circular Plates 


The differential equation for the free transverse vibrations of a 
circular plate is given as 



i 

r br 


-f- 




2 

w 


+ p!1 ^ . 0 


(5.23) 


where 


D 

P 

E 

Ut 

8c 

h 

w(r^ 0, t) 


Eh' 




flexural rigidity of the material 


mass density of the material 
Young's modulus 
Poisson's ratio 
gravitational constant 
plate thickness 

displacement of a point on the middle surface of 
the plate„ 


Equation (5n23) is based generally on the assumption of equation 


(5.2). 
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(a) 


(b) 


Figure 5.2 Clamped circular plate and k-space 

The method of separation of variables can be used to solve 
equation (5 „ 23) c 

Assuming that the solution of equation (5,23) is 

w(r, 0, t) = W(r, 0) (5 

and substituting equation (5.24) into (5.23)^ yields 

1 Bw ,1 B^W,^ 


, 2 , 



where 


.4 _ (JO 

A ~ r) 


(i^ - a constant 




12(l-p. )p 


Again applying the separation of variables method for solving 
equation (5.25) gives 


W(r,e) = R(r) ^0) . 


(5.26) 


Substituting equation (5.26) into (5.25) gives 



- ^ - (±\^ 
r dr 


2 

+ — ^) R = 0 


(5.27) 


d 6 1 

^ + n^^t =0 (5.28) 

d0 

2 

where n is a constant. Solving the equations (5.27) and (5.28) 
and substituting into (5.26) yields 


W(r,0) = [A^J^(Ar) +B^J^(iXr) +C^Y^(Ar) +D^Y^(iXr)] 

[E cosn9+F sin n6] n = 1. 2. 3, ... (5.29) 

where A , B . C . D , E and F are the arbitrary constants and depend 
n'' n'’ n-^ n-^ n n ^ 

on the boundary conditions of the plate. 

J and Y are the Bessel functions of the first kind of order n 
n n 

and the Bessel functions of second kind of order n respectively. 
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Since the functions Y (r) become infinite for r = 0 

and D must be zero in the case of a solid plate, 
n ^ 


(5,29) becomes 


, the constants C 
^ n 

Hence equation 


W(r^0) = [A J (Xr) + B J (iXr)][E cosn0+E sinn6] . (5.30) 

^ ^ ^ n n n n' n n ■' 

For the plate under consideration, the boundary conditions are 

W(r,0) = = 0 at r = a . (5.31) 

By applying the boundary conditions to equation (5.30), the frequency 
equation for a circular plate can be expressed as 


X 


m, n 



(5.32) 


where uj is a natural frequency corresponding to the mode charac- 
terized by the eigenvalue X and 8 is a constant defined as 

m, n 


Eh^g, 


P = 


)p 


Thus, the exact angular frequency is given as 


u) = X' 


2 

Eh^g 


nij ■ 1 o / 1 

^ I2(l-|i ) p 


rad/sec 


(5.33) 


Now, for large values of X, 


(5.34) 
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Therefore^ 
plate with free 


for high frequencies^ the frequency equation for a 
edge or clamped edge conditions is given as 


U) = — ^ (m + 2n) y/ 

43“^ Vl2(l-n)p 


(5.35) 


KC (m + 2n) ^ 

/ ^ Li 

4a 


m, n = ly 2, 3;, 


(5.36) 


where 


a = radius of the plate 
K = radius of gyration 


ion 


12(U/) Vl2 


C = longitudinal velocity of wave propagation. 

Lj 


The frequency equation can also be written as 


(JD 


. mn 2 

4a a ^ 


KC. 


(5.37) 


Defining the wave numbers and as 


, miT j , nTT 

4 = 2 l *^2 = T ^ 


(5.38) 


the change in the wave numbers from one mode of vibration to the next 
is given as 


Aki = ^ and 



(5o39) 


The wave propagation, in the case of a circular plate, takes place 
In the two directions so that the k-space is also two dimensional. The 
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equation for the number of resonant frequencies is given as 


N(u>) 


Ak 


J J* 


(5.40) 


In this case, cylindrical co -ordina tes can be utilized to inte- 
graiie over the surface of k-space. Letting k^ ~ r cos 0 and k^ = r sin 0 
and substituting the values for the change in wave numbers and convert- 
ing CO cylindrical co-ordinates, equation (5.40) becomes 


N(od) 


2 


-r I 


(5.41) 


The frequency equation (5., 37) can be written as 


Changing 

2 . 

r gives 


(kj . k2 

equation 



(5 42) to 


cy i indrica 1 


( 5 „ 42 ) 


co-ordinat es and solving for 


2 

r 


(JD 1 

K C ^ ~ 2 

L ( sin 0 + cos 0) 


(5..43) 


2 

Substituting the value of r into equation (5„41) and carrying out 
integration over values of 0 in the quadrant 0 < 0 < ^ for which the 
integrand is real, gives 


N((«) 


2 

a oj 


TT L 


(5.44) 
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U)3 

Defining v = , equation (5.44) reduces to 


av 

N(V) = - 5 - 
n K 


(5.45) 


where v is a dimensionless frequency. 

Differentiating equation (5.45) with respect to v gives 

n(v) = . (5.46) 

TT K 


In case of a circular plate^ there are two modes of vibration 
for each frequency^ hence the modal density of a circular plate is 
doubled for any frequency and expression (5.46) reduces to 

n(v) = . (5.47) 
TT K 


This is an expression for the modal density of circular plates. 

If the plate thickness is h^ radius of gyration is h/\fl% then 
expression is given as 


n(v) 


4a\[3 


(5.48) 


5o4 Discussion 

From the expressions developed for the modal density of flat 
rectangular and circular plates^ it can be concluded that the modal 
density of a flat plate is a constant for a given plate and thus is 
independent of frequency. 
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Now for a rectangular plate for each frequency^ there is just one 
mode of vibration^ whereas for the circular plate^ there are two modes 
of vibration for each frequency. Therefore for two plates of equal 
area^ thickness and of the same material, one circular and the other 
rectangular, the ratio of modal densities is found to be 


n(tu) j. 


TT hCj^ 


2h TTC 

Li 


8 

"7 

n 


1 . 


Thus for a given frequency (high frequency due to an assumption in 
the circular plate derivation), the modal density of the rectangular 
plate will be approximately equal to that of the circular plate. 
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6, MODAL DENSITY OF THIN CIRCULAR CYLINDERS 


6o 1 Introduction 

In this chapter^ thin walled circular cylinders are considered and 
the expressions for the modal density developed in three different ways 
are discussed in detail. The first presentation is that of Heckl 
(1962) in which the expressions are found for the cylindrical shell 
alone. The second representation is that of Bolotin (1963) in which 
the general shell of evolution is discussed and then applied to the 
case of the thin cylindrical shell. The third representation is 
essentially a modification of Bolotin's work for the specific case of 
the cylindrical shell. 

Shells simply supported at their edges are considered for 
developing the expressions. However the effects of boundary conditions 
on the vibrational modes are limited^ and hence the edge conditions 
are of little significance in the modal density expressions except 
for the first few modes. 

All the three representations are discussed in the literature 
(Miller and Hart_, 1967) and they are reproduced here and expressed as 
a function of dimensionless frequency. The graphs of normalized model 
density versus frequency are plotted. 

6,2 First Representation 

For a thin infinitely long cylindrical shelly the equations are 
given as 

2 

cvV + n + LlkaV =iVPg/(JDph (6.1) 

o t a o c' 
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( 6 . 2 ) 


n^C + [i1q + +%(l-|j,)k.^a]Vt +%(l+|J,)nokaVa = 0 

p,kaV +%(l+|i)nQkaVj. +[k^a^ +%(l-|j,)n^ - v^]V^ =0 (6.3) 

where 

My V and V are all radial^ tangential and axial components 
t a 

of velocity amplitude 

n^ = the half number of modes in the circumferential 
direction 

k ~ wave number in the axial direction 
\i - Poisson's ratio 
h = shell thickness 

- gravitational constant 

0 ) = frequency of vibration 

- amplitude of excitation 
a = cylinder radius 

V - - dimension frequency of vibration 

a = velocity of wave propagation in the shell material 
a = 1- +{(n^ +k^a^)^-^[np(4...^.) - 2 *^i](l..-^)'^)hVl2a^ . 

Equations (6.1), (6.2) and (6,3) can be solved for the impedance 

of the cylindrical shell and letting the impedance go to zero^ the 
following frequency equation can be obtained- 


2 

V 


(l-v^) (!!^)^x [(! 1 ^)^ + n 2]-2 


%[nQ(4-n) - 2 - n](l-n)-h — 


12a 


(6.4) 
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The approximate frequency equation obtained by neglecting some of 


the terms is given as 

2 r 2 , 2 ^ 2,-1 ^ 2 ^ 2 , 2-,2 

V = [a (Hq +cr ) + P(np +0- ) ] 


( 6 . 5 ) 


where 


O' = 


mrra 


P 


2V3; 


The terms neglected have little effect on the frequency expression 
for frequencies above the ring frequency (v = 1) » However^ below the 
ring frequency the effect may be as much as forty percent of the actual 
value. 

Now solving for a and then summing over all possible values of n, 
the equations as obtained by Heckl for the number of resonant frequen- 
cies and the modal density are given as 


n 

N((v) = S ^ 


( 6 , 6 ) 


and 


n(w) = 


rra 


J 




0,1 


ba 


dv 


( 6 . 7 ) 


where the lower limit is 1 for v < 1^ and 0 for v > 1. 
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Simplifying further results into the followings: 


For V > 1 


N(u)) 


V3 t ao) 

2C_h 

J-i 


n(u)) 


^3 I a 
-2^ 


For V < 1 

N(to) = [%(2v-l) [%n + arc sin (2v-l) + (v.^v^) 
n(o)) = [%TT + arc sin (2v-l) ] 


( 6 . 8 ) 


(6.9) 


( 6 . 10 ) 

( 6 . 11 ) 


These are the final expressions obtained by Heckl for the number 
of resonant modes and the modal density of a thin cylindrical shell. 

These expressions as functions of dimensionless frequency are 
given as: 

For V > 1 


N(V) 

\[3 tv 
2h 

(6.12) 

n(v) 

][3l 

2h 

(6.13) 


For V < 1 


N(v) 


3lv^^^ ^ 3yj I 
Strap 4Trah 


n(v) 


9V3-tv^/^ 

Sirah 


(6.14) 
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6.3 Second Representation 


The diff. equation of motion for thin shell of revolution is 
given as 


DA Aw _ 







Eh 


AA^ 


Ro 



d w 



0 


(6.15) 


(6.16) 


where 

and are the general curvilinear co-ordinates 
and are the principal radii of curvature 
3 3 

T. Eh Eh 1 ^ 

D = — » jy- = the plate stiffness 

12(l-n^) 

p = density of the material 
h = thickness of the shell 
w = normal deflection 

^ = stress function for the middle surface 

E = Young's modulus of elasticity 
O) = frequency of vibration 
g^ = gravitational constant. 


The solution of the equations (6. 15) and (6. 16) gives the follow- 
ing frequency equation: 


2 

U) = 


Eg 


c r /T 2 ^ , 2. Eh 


(kjX+k^) 

2 2~ 

(kj . k2) 


(6.17) 


40 



where 


where 

Now 


"1 

the 



and 

2 

ITITT 


and 

number 


are the wave numbers and are given as 

and = — m,n =1^ 2^ 3^ o,. (6.18) 

Z 32 

are the principal dimensions of the shell surface, 
of resonant modes in the shell is given as 


= Atdiac -T J" • 


(6.19) 


12s 


The change in the wave numbers Ak^ and Ak^ from one mode of 
vibrations to the next is given by 


Ak, = — and Ak^ = — (6.20) 

1 3i 2 32 

where and a^ are the principal dimensions of the shell surface^ 

Substituting the values for the change in wave numbers and con- 
verting to cylindrical co-ordinates^ equation (6.19) reduces to 


12 

N(U)) = 

J J rdrde 


n 

S 


= ^1^2 

e„ „ 

J ' ■ 

(6.21) 

~2t? 
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(a) 


(b) 


Figure 6,1 Generalized rectangular region and k-space 


Converting the frequency equation (6. 17) to cylindrical co- 

2 

ordinates and solving for r , gives 


r = [cu -n^(Xcos e + sin 6) J- (gg-) 


( 6 . 22 ) 


where 


1 ^ 


Thus the expression for the number of resonant modes becomes 


a a 

_ 1 2 ph r r 2 ^ 2 , 


N(U}) = — L £ J [oj"^ - Q^(Xcos^6 +sin^0)^] ^ 


2tt‘' ^^c “©^(u)) 


de 


(6.23) 


42 



Equation (6.23) is an integral expression for the number of 
resonant frequencies up to a boundary frequency for a thin shell of 
revolution. The integration is taken over the values of 6 in the 
quadrant 0 < 0 < for which the integrand is real and positive. 

Differentiating equation (6.23) with respect to frequency U) under 
the integral sign using Leibnitz’s rule gives 


n(o>) 


^1^2 


2tt 



d0 

[cu^ -n|(Xcos^ e +sin^ 6)^] ^ 


(6.24) 


This is an expression for the modal density of thin shells of 
revolution. 

These equations as written by Bolotin can be represented in the 
form as shown: 


N(u)) 


,ph .h u A 

(^ — ) (W H (— 

4rr Dg v 

c 


X) 


(6.25) 


n(uj) 


!i!i H fi X) 

4tt ■’ ’ 


(6.26) 


where 



V CD tjuR^ 

62(00) 

H(i , X) = i J [1 - i (Xcos^ 0 +sin^6)^]^ de 

6j^(a)) \o 




x> -ll 


02 


d0 


0^(00) 


V 


( X cos^ 6 + sin^G) 
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Now for cylindrical shell in particular 
X = 0 

and a^ are the dimensions of the shell surface 

hence = t - length of cylinder 

a2 = na - one-half the circumference of the cylinder. 

The reason for taking only half of the cylinder into account is 
that the cylinder is a closed surface and that the vibrational modes 
are limited to one-half by that fact. 

Hence equations (6,25) and (6.26) take the form 

N(co) 

°C 

n(c«) = ^ 0) . (6.28) 

Now rewriting the expression (6.27) as a function of dimension- 
less frequency v gives 

\Tr i n 

N(v) = H(- , 0) (6,29) 

and 

\/T 1 

n(v) = -^ 0) . (6.30) 

The expressions (6,28) and (6.30) are expressions for the modal 
density of a cylindrical shell of length t and radius a. 
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The expressions for may be used in the elliptical integral form 
as follows: 


For V > 1 


Hi(- , 0) 






(6.31) 


For V < 1 


Hi(i 0) = — (\ T ) 


rr\/T/v 


(6.32) 


where represents the complete integral of the first kind. 

6.4 Third Representation 

As stated in Section (6.4), the number of natural frequencies for 
a thin shell of revolution is expressed as 

0 „ 


N(tu) = J [ou^ -n^(Xcos^ e +sin^ 0 )^]'® de . (6.33) 




2tt" ■ 0 


However for a thin cylindrical shell in particular 
a^ = -L = length of cylinder 


a^ = an = half the circumference of cylinder 


X = 0 


Then equation (6.33) reduces to 

e. 

, 


N(0)) = ^ (P|-)^ ; "[u,2.0^sinS]^ d0 


e. 


(6.34) 
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where the upper and lower limits of the integral are taken in the first 
quadrant (0 < 0 < j) in such a way so as to keep the integrand real and 
positive. 

Rewriting equation (6. 34) as a function of dimensionless frequency, 
it reduces to 


N(V) 


t\f3 

nh 



[ V - sin 6] 


d0 


(6.35) 


The upper limit on the integration holds for v < 1. For v greater than 

TT 

or equal to one^ the upper limit ^ is used. 

Differentiating equation (6.35) with respect to v gives 


n(uj) 


2h 


- J 

TT ‘J, 


sin'^y^ 


[1 


de 

sin^e]^ 


(6.36) 


Again the upper limit must be ~ when v is equal to or greater than 

one. 

Equation (6.36) is the expression for modal density of thin 
walled circular cylinders and is referred to as the modified Bolotin’s 
result. It can be evaluated numerically by means of Simpson's rule 
using a digital computer. 


6.5 Graphical Results and Discussion 
Figure 6 . 2 shows the variation of the modal density for the 
number of natural frequencies for three different representations. 
Above the ring frequency (v > 1)^ all three representations give the 
identical results as v becomes very large. However below the ring 
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frequency the differences are quite noticeable. The results of Bolotin 
and the modified Bolotin analysis are slightly different for the modal 
density. The reason for this slight difference in the modal density 
curves is due to an approximation made by Bolotin in order to express 
the modal density in elliptical integral form„ The results of Heckl 
below the ring frequency are lower than other results^ but since it 
was shown earlier that Heckl ‘s results would be on the conservative 
side^ it is reasonable to assume that this is the reason for the 
difference. Comparison of expressions for modal densities with that 
of a plate shows that modal density of a thin cylindrical shell above 
the ring frequency (v > 1) is equal to one-half the modal density of a 
flat plate with the same surface area. 


48 



7. MODAL DENSITY OF THIN SPHERICAL SHELLS 


7 . 1 Introdviction 

In this chapter, the thin walled spherical shell is discussed and 
the expression for modal density is developed in two different ways 
following the second and third representation used in deriving the 
expression for model density of thin circular cylinders. In the first 
representation_, the general shell of revolution is considered and then 
applied to the case of the thin spherical shell, whereas in the second 
representation the spherical shell is considered in particular.. 

The frequency equation derived for the thin shell of revolution is 
in general for the shell with simply supported edges; however, the 
expressions for the number of modes and modal density hold good for 
all the boundary conditions, since it is examined in detail that for a 
spherical shell, edge effects never dominate the mode shapes (Bolotin, 
1960) . 

The expression for modal density is obtained in terms of dimension 
less frequency and a graph of dimensionless medal density versus 
frequency is plotted^ 


7o2 First Representation 

As stated in the previous chapter, the expression for the number 
of resonant modes for free transverse vibrations of a thin shell of 
revolution is given as: 


N(cd) 


^^2 


2n 



n^(xcos^e-rsiti^e)^ de 


(7.1) 
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where 


and a^ are the general curvilinear co-ordinates 
and are the principal radii of curvature 
R. 


X = 


Rn 


p = density of material 
h = thickness of the shell 
E = Young* s modulus 


g^ = gravitational constant 

(JD = natural frequency of vibration 


1 ^ 
^ p ^ 


Differentiating expression (7.1) with respect to frequency cu under 
the integral sign using Leibnitz's rule gives 


a, a 


0. 


n(cju) = 


1^2 ,ph % p " de 

2^2 ^Dg^' ^ 2 , 2 ^^_.^ 2 ^, 2 ,% 


(7.2) 


0^ [a)"^O^(Xcos"0+sin"0)^]^ 


This is an expression for the modal density of a thin shell of 
revolution. These expressions (7.1) and (7.2) as written by Bolotin 
can be represented in the form: 


N(u)) = 


a a 
12 

4tt 


,,x ) 

c 


(7.3) 


n(uj) 


^1^2 ph k 

4tt ^Dg ^ 

c 



(7.4) 


50 


where 


H(i , X) = - r [1 - i ( Xcos^e + sin^ 6 )^]^ de 

^ e^((0) V 


H 


d0 


■’ IT i*Q^ [1 _ -^ ( Xcos^e H-sin^e)^]"^ 


V 



dimensionless frequency . 


For a spherical shell of radius a 


X = 1 


a 


1 


rra = half the circumference of sphere , 


The reason for taking only half of a sphere into account is that 
the sphere is a closed surface and that the vibrational modes are 
limited to one-half by that fact., 

Hence equations (7„3) and (7«4) take the form 


N(cu) = 


2 

rra 

“2T 




1) 


(7,5) 


and 


n(oo) 


2 

Tia 

~7T 


/ Ph - v% 


H,(i 

l^V 


1 ) - 


(7.6) 


Again, rewriting the expressions as a function of the dimension- 
less frequency v gives 
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N(v) 


(7.7) 


= I|V3 V H(1 , 1) 
2h 


and 


n(v) 


~2h~ 



1 ) 


(7.8) 


Now 


Hi(- . X) 




n 


de 


0^ [1 - ( X cos^0 + sin^ 0) 

V 


(7-9) 


The integration with respect to 0 is carried out over that part of 
TT 

the quadrant 0 < 0 < -^ which the integrand is positive and real„ 

For a spherical shell X = 1 and the integral (7=9) can be expressed as 


Hi<i , 1) = 0 

(V < 1) 

(7.10) 

"i4 ' “ — ~ 

i V f-s 

(v > 1) . 

(7.11) 




Hence equation (7-8) is written as 

n(v)=^^-^ (V > 1) (7.12) 

n(v) =0 (V < 1) . (7.13) 

Equations (7.12) and (7.13) represent the modal density for a thin 
walled spherical shell. 
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7o3 Second Representation 


The expression for modal density can also be obtained by finding 
the number of resonant frequencies for a spherical shell and then 
differentiating with respect to frequency^ in the same way as the 
modified Bolotin results for cylindrical shells. From equation (7.1)^ 
the number of resonant frequencies for a thin shell of revolution is 


a a 

N(uj) = I J [u)^ - n?( Xcos^e + sin^9)^]^ d6 . 

2n ‘^e^(ou) ^ 

(7.14) 


Now for the thin walled spherical shell 


a^ = a^ = na = half the circumference 


X = :^ = 1 . 
^2 


Hence expression (7.14) reduces to 

^2 2 

r,., 


N(UJ) = ^ J ■ 

c 0- 


(7.15) 


Converting equation (7,15) in terms of the dimensionless frequency 


gives 


N(v) 


iEj vCl-4]^d6 


(7.16) 


Here again the limits on the integral are taken in the quadrant 


0 ^ 0 < 5 such that integrans is real and positive. 
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Hence the following integral expression for the number of resonant 


frequencies is obtained: 
N(v, . ip / 


sin 


_ 1]^ de . 


(7.17) 


The upper limit on the integrations holds for v < 1- For v > 1 
the upper limit ~ is usedo 

Differentiating with respect to v 


n(v) = ^ 


J 


sin"^V^ 




de 


(7.18) 


Now for V < 1 the integrand is negative^ hence 


n(v) = 0 


(V < 1) 


(7.19) 


and for v > 1^ using the upper limit — , it gives 


n(v) 


V^rra v 


(V > 1) . 


(7.20) 


Equations (7^19) and (7.20) represent the modal density for a 
spherical shell. 


7o4 Graphical Results and Discussion 
Expressions obtained for modal density following two different 
approaches gives identical results. In plotting the graph_, the expres- 
sion is normalized so as to make it independent of geometry. 

Figure 7.1 shows that modal density has a singularity at v = 1, 
below which the modal density is zero. For v > 1^ (above the ring 
frequency) the modal density of the shell decreases monotonically and 
as37mptotically approaches that of the flat plate as v becomes very large. 
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Figure 7.1 Normalized modal density versus dimensionless frequency for 
thin spherical shells 
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8. MODAL DENSITY OF THIN CONICAL SHELLS 


8,1 Introd uction 

In this chapter^ the thin walled conical shell is discussed and 
the expressions are developed for the number of resonant modes and 
modal density. These expressions are applicable to a wide range of 
geometry and are valid over a frequency range sufficiently wide to be 
of engineering value. 

The expressions for the number of resonant modes are developed 
from two separate frequency equations for thin conical shells using the 
k-space integration technique for frequencies below the lower ring 
frequency and a different approach for frequencies above the upper 
ring frequency. The expressions obtained for modal density for two 
cases are normalized and plotted versus the dimensionless frequency. 

The problem of conical shells was readily available in the 
literature and is reproduced for graphical representation (Miller, 

1969), 


8,2 Frequency Equation One 

The frequency equation for a thin conical shell is given as 


4 

D r ^ 


2pC,ELh SLn\k n 

(U = [__] 1 


sin t|; n 


+ ( 


^ / 2 

m 4m 


( 8 . 1 ) 


.4, .2, 

sin \|f sin \j; 


tan t|; 


8a 


4 2 ---- - 

+ ( 7 . o^) ^ (l-cy^) ^ (l-cv^) + — ^ (1-cy^) ] 


sin^\jr sin^ijf ^ 


2a 


n 


4a 


n 
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where 


U) 

E 

P 

D 


M. 

h 



L 

m 

n 

a 

n 


natural frequency of vibration 
modulus of elasticity 
gravitational constant 
density of shell material 


stiffness of shell material 


Eh'^ 


12 ( 1-0 


Poisson* s ratio 

the thickness of shell wall 

one-half cone angle at apex 

. 

truncatLon ratio = 

i-r 

length of cone truncation, apex to top slant length 

length of cone^ apex to base slant length 

number of circumferential waves 

number of one-half longitudinal waves 
nn 


In deriving the above frequency equation, the following assump- 
tions are made: 

(1) The circumf erential modes are independent of the longitudinal modes- 

(2) The modes in the circumferential direction are sinusoidal and 
uniform over the length of the cone- 

Equation (8-1) may be written in dimensionless form by defining a 
dimensionless frequency ratio and longitudinal and circumferential wave 
numbers in the following manners 


2 2 

2 00 pL 

\) = z 


k_ = a 
1 n 


m 

sin^ 


( 8 . 2 ) 
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Substituting equations (8.2) into equation (8.2) and rearranging 
terms results in the following expression for the dimensionless 
frequency: 


2 2 
12L (1-V ) 


/ / 9 9 


l-“i / 9 2 


4 




l-.Q', 


10 6 


-) + (k 


4k^) 


1-cv^ 


1-cT 

2 4 1 

+ cot \|;[k^ — g— 


2 3(l-aJ) 

—8 


] 


X [k" 


1-Q^ 

1 

10 


- k 


2^-i 


3(1-“ ) 
+ — 


(8.3) 


Now the number of resonant modes for a thin conical shell is given 
by the double integral 

■ 4k^ S J' W-*' 

where the integral is to be taken over that portion of the first 
quadrant where the k-space exists^ 

The integration of the equation (8.4) to obtain an expression for 
the cumulative number of eigenvalues or the resonants modes in the 
usual way is a somewhat impractical approach to the problem and hence 
the number of eigenvalues above some selected frequency is obtained 
utilizing a different approach* 

It will be of some value to first define upper and lower ring 
frequencies. The upper ring frequency is defined as the frequency at 
which the longitudinal wave length is equal to the circumference of the 
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small end of the cone. In dimensionless form the upper ring frequency 
would be given by 


V (upper ring) 


1 

a^sinijf 


(8.5) 


This is equivalent to (Jo times the small radius of the cone divided by 
the longitudinal wave velocity equal to unity. 

Similarly the lower ring frequency is defined as the frequency 
at which the longitudinal wave length is equal to the circumference 
of the large end of the cone. In dimensionless form it is 


V (lower ring) - — : — - 
sinilf 


( 8 . 6 ) 


This is equivalent to U3 times the large radius of the cone divided by 
the longitudinal wave velocity equal to unity. 

The frequency equation (8.3) shows that dimensionless frequency is 
affected by these geometric parameters assuming that Poisson's ratio 
is constant and equal to 0.3. These are the cone angle (ilr)^ the 
thickness over length ratio (^) and truncation ratio (cv, ) . 

Li i 

Using frequency equation (8.3)^ the eigenvalues may be computed 
for different values of m and n^ the circumferential and longitudinal 
wave numbers respectively. In this manner for different cone 
geometries^ the number of eigenvalues occurring above certain specified 
dimensionless frequencies may be obtained using digital computer and 
results can be plotted in a graphical form. 

Now the results obtained by normalizing with respect to cone 
geometry in the frequency range above the upper ring frequency (Miller, 


59 



1969) can be expressed in the form 




(8.7) 


The graphs of F(v) versus dimensionless frequency can be plotted 
varying various geometry parameters and it can be found that the value 
of F(v) is independent of the cone geometry above the upper ring 
frequency (Miller^ 1969) and is given as 


F(v) = 2.0v 


for 


V > 


1 

cy^sin^ 


( 8 . 8 ) 


The number of resonant modes is given by 


N(v) 


2.0 


L sin i|f 

[ 


rfh 



V 


(8.9) 


Differentiating equation (8.9) with respect to v gives 


L sin \|f 

n(.) = 2.0 [ ^ 


4/5 


-] 


( 8 , 10 ) 


This is the expression for modal density for thin conical shells 
above the upper ring frequency. 

8.3 Frequency Equation Two 

The second frequency equation for a thin conical shell is given as 


^ 4 l.a, 2 1-a^ 3(l-aJ^ 2 

D . m 1. n r 1 ^ 1 T 

2, ^.4,^ 2^ 2. ^"“8 7~2 ^ 


2 , 8 

tan i|r 


8a, 


^ g E EL h sin T|f 

U) = [-^] [ J ^ 

pL 4 1-c.. l-< l-a, 3(l-a.) 


] • 


( 


\) ( i i + ^ ) 

' ^ -in n A ^ 


4 210" 2 4 

sin^ilf 2a 4a 

^ n n 


( 8 . 11 ) 
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II 


The notations used in the above equation are the same as those for 
equation (8ol)o In deriving the frequency equation (8.1) it is 
assumed that: 

1. Mode shapes are axially symmetric^ and 

2o Longitudinal bending is small when compared with circumferential 
bendingo 

Defining the dimensionless frequency parameter^, the circumferen- 
tial wave number and the longitudinal wave number in exactly the same 
way as in Section 8.2, equation (8.11), in the dimensionless form can be 
written as 


2 

V 


327(^2—) rr-J 

EL h 8k, 


l.„ 3a.c^> 

^2 2 ^ 10 2 ? ^ 

2k^ 


( 8 . 12 ) 


Now the number of resonant modes or eigenvalues for thin conical 
shells is expressed in the double integral fom as 


N(v) 




(8.13) 


where the integral is to be taken over that portion of the first 

quadrant where the k-space exists. 

The region over which the integral (8.13) is to be evaluated is 

bounded by an upper and lower value of k^^ referred to as b and a 

respectively. The region is also bounded by upper and lower curves 

which will be referred to as [k«] and [k^]« respectively. Equation 

Z u Z i/ 

(8.13) therefore can be written in the following form: 
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»<''>■ 4i^ l" '“‘'i «•“> 

1 Z a 

where a and b are functions of v and the cone geometry parameters and 

[k ] and [k ]. are functions of v, k- and cone geometry parameters^ 

Z U Z t/ X 

In equation (8c 14) Ak^ and Ak^ are the changes in the longitudinal 
and circumferential wave numbers respectively from one mode to the next 
and are given as 


Aki = — ^ 
1 1-Q', 


Ak_ = — : — - 
2 sini|f 


(8.15) 


Substituting equation (8.15) into (8.14) and rearranging the 
terms gives 


1 a 






dk. 


(8.16) 


The upper and the lower bounds of the space [k ] and [k^]» may 

Z U Z i/ 

be obtained from the solution of the equation derived from frequency 
equation (8.12) and may be expressed in the following form; 


and 



(8.17) 


(8.18) 
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where 


A = 


B = 


C = 


D 

EhL^ 



v^[-16k^(l +o^ +a-j) +8k^(l +a^ +0-^) - 120] 

5 cot^ilt[k^(l + a^ +0-^ +Qf^) - 3k^(l +0'^^)]^ . 


Therefore equation (8. 16) may be expressed in the following form: 


N(v)[ 



1 ‘ 


(8.19) 


The upper and lower limits of the space, b and a, can be obtained 
by equating the upper and lower bounds of the space given by equation 
(8.17) and (8,18). The expression resulting from this process is given 
as follows: 


SkJ + 


Tk, + U = 0 


( 8 . 20 ) 


where 

S = [-3(l-a^)/(l-aJ)] - [4v(l-a^)L 5h cot tCl-aJ) ] 

T = 4v(l-o^) L ^3(l-|x^) / h cot (|;(1-Q'^) 

U = -6v(l-ap L ^3(1-11^) / h cot iKl-o^) . 
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For the truncated cone^ three real roots are obtained; the largest 
root is denoted as b and the second largest root is denoted by a. For 
the closed cone^ one real and two imaginary roots are obtained. The 
real root is denoted as b and the real part of either of the imaginary 
roots is denoted as a. 

The evaluation of the equation (8.19) for the number of eigen- 
values and also upper and lower limits of the k-space is handled by 
numerical procedures on an IBM 360 model 75 digital computer and graphs 
are plotted. 

Now from the graphical results (Miller, 1969) it can be concluded 

that the number of modes (N — : — — ^ ^) varies with changes in cone 

sin^(l-Ck'^) 

V h 

angle directly as (tan;|;) with changes in thickness inversely as — ^ 

1/4 

and with changes in truncation inversely as (l-cv^) ^ . Hence the 

number of modes may be normalized in the following manner: 


N(v) 


n 

sin\jf( 1-Qf^) 


L(tani|f) 


1/4 

Tfl ^ 


G(v) . 


( 8 . 21 ) 


In deriving frequency equation two^ it was assumed that the 
contribution to the differential equation due to longitudinal bending 
is small in comparison with the contribution due to circumferential 
bending which limits equation two to the lower frequencies. Hence, 
the results of the analysis based on frequency equation two are valid 
below the lower ring frequency. 

The graphical representation reported by Miller (1969) indicates 
that G(v) is independent of the geometry of the conical shell except the 
variations in the vicinity of values associated with N(v) = 1 and is a 
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function of only the dimensionless frequency v. Moreover_, the 
normalized number of eigenvalue curves are straight lines on log-log 
paper except near the values associated with the first few resonant 
frequencies. Hence equation (8.21) may be approximated as 


N(v) 


TT 

sinil/Ci-cv^) 


[ 


L(tani[f) 


1/4 

1/2 ^ 


0,876V 


3/2 


( 8 . 22 ) 


This equation is valid for the frequencies below the lower ring 
frequency only. 

Hence the number of resonant modes for a conical shell below the 
lower ring frequency is given as 


N(v) 


0.876 


[ 

sin\|; 


X 


h(l ^ ^2 


(8.23) 


Differentiating equation (8.23) with respect to v gives 


n(v) 


1.31 [ 


TT 

sin\|; 


X 


1/2 

172 ^ 

L(tan)|f) ' 


(8.24) 


This is an expression for modal density of thin conical shells below 
the lower ring frequency. 


8.4 Graphical Results and Discussion 
Expressions (8.10) and (8.24) represent the modal density for a 
thin conical shell and variation in modal density is shown in Figure 

8 . 1 . 

In plotting the graphs the expressions are normalized so as to 
make them independent of geometry. For the variations below the lower 
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ring frequency cone angle appears as a parameter as the lower ring 
frequency depends on the cone angle. The upper ring frequency has not 
been indicated in the figure since it is a function of the truncation 
ratio of the cone and may vary anywhere from the lower ring frequency 
for a completely truncated cone to infinity for a closed cone. 

The graph and expressions are invalid for cones with large cone 
angles as well as with little or no truncation at all as the shell in 
these goes into a so called plate mode (above upper ring frequency 
solution) . However it can be concluded that expressions are valid over 
a wide range of cone geometries and frequency ranges of practical 
interest. 
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9o MODAL DENSITY OF COMPOSITE STRUCTURES 


9c 1 Introduction 

In this chapter the additive property of modal density for 
composite structures is discussedc Modal density of certain basic 
structures such as rods^ beams, plates^ cylinders and spheres have been 
already discussed in the previous chapters. However these basic 
structures rarely occur in a real application in engineering as 
separate elements. Therefore the modal density of composite structures 
must be considered. 

The composite structure analyzed consists of two beams joined at 
right angles to form an L-shaped frame. The case of composite struc- 
ture was readily available in literature (Hart and Desai^ 1967). 


9c 2 Composite Structures 

A composite structure is composed of a number of substructures 
which may be taken as the sum of the basic structural elements. 
Assuming that the modal density of each substructure is known^ it is 
postulated that the modal density of a composite structure is equal to 
the sum of the modal densities of its components. 

If the component of the composite structure exhibits modes 
within the frequency interval Acu^ then its modal density at the center 
of the band Au) is defined as 


n^(CA)) 



(9.1) 


Thus. the modal density of the composite at O) would be given 
by 
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m 


(9.2) 


where the sununation extends over the total number of elements^ m, that 
give rise to the composite structure. 

To demonstrate analytically^ a composite structure consisting of two 
beams joined at right angles to form an L~shape is considered (Hart and 
Desai^ 1967). The subsystems may then be supposed to be two beams as 
illustrated in Figure 9.2. 


Figure 9.1 Composite structure 





qi qi 

Figure 9.2 Two substructures 
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It is assumed that = 't and that both members are of the 
same material. 

If the additive postulate for modal density holds good^ then the 
total numbers of resonant modes for the composite structure are given 
as follows: 


Na(U) = Ng(X^) + N^(U) . 


(9.3) 


By considering the frequency equation for the composite structure^ 
a graph of against X-t can be plotted. The relationship derived 
from this graph (Hart and Desai^ 1967) can be expressed as follows: 


N^CA.'C.) = ^ (U) 


(9.4) 


Now the resonant frequency for a composite structure derived from 
the frequency equation is given by the expression 

2 


OD 


(Xl)‘ 


(9.5) 


where 


= longitudinal wave velocity 
K = radius of gyration. 

Equation (9.5) gives 


(U) 


' P T/- 




(9.6) 


Defining dimensionless frequency v 
u>t 


V = 


(9.7) 
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Equation (9.6) can be written as 


■ 

Substituting (9.8) into (9.4) gives 


(9.8) 


N^(v) 


16 

25 '' K 


(9.9) 


Differentiating equation (9.9) with respect to dimensionless 
frequency V gives 


Ha(v) 


_8 iTTJ 

25 ’ Kv 


(9.10) 


Equation (9.10) gives the expression for modal density of a 
composite structure. 

The modal density for a beam^ irrespective of the boundary condi- 
tions^, can be written as 


n(v) 



(9.11) 


Since the composite structure is constructed of two identical 
beams^ the sum of the modal densities of the substructures is 


,^(V) . »c(v) . i . <9.12) 

Comparison of equations (9.11) and (9.12) proves that additive 
property of modal density holds good for this particular composite 
system. 
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9.3 Graphical Results and Discussion 


The additive property of modal densities for composite structures 
holds for the composite structure composed of two identical beams 
welded at right angles as verified analytically. 

Figure 9.3 shows that modal density of a composite structure 
varies along a straight line on a log-log scale and variation with 
respect to frequency is proportional to the beam except that the 
magnitude of modal density is doubled. 
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10 . MODAL DENSITY OF SHALLOW STRUCTURAL ELEMENTS 


10 . 1 Introduction 

In this chapter the expressions for the modal densities of some 
shallow sandwich shells^ orthrotropic plates, pretwisted plates^ 
plates subjected to in~plane forces and shells on an elastic foundation 
are presented. 

The expressions are developed on the basis of shallow-shell theory 
and neglecting the effect of longitudinal inertia and hence only the 
frequencies of transverse vibrations are considered. The coupled 
longitudinal modes cannot be obtained from these expressions. However 
the effect of the longitudinal modes on the modal density of a shallow 
element is negligible as these modes occur only at widely spaced 
intervals over the frequency spectrum. All the expressions are 
strictly valid only for simply supported structural elements, however 
for large values of U) it is reasonable to suppose that asymptotic 
relations for modal density are relatively independent of the boundary 
conditions (Courant, 1953) . 

All the structural elements considered are discussed in detail 
(Wilkinson, 1967) and they are reproduced here for the graphical 
representation. The graphs for modal density versus frequency are 
plotted for various elements for specified dimensions and properties. 

10.2 Sandwich Shells 

The governing equation of motion of the shallow sandwich element 
of constant curvature is given as follows; 
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"y . 


n ] w = 0 


4Ti^(Pj^h^ +P2h2)i»^^ 


2L.2L 

2 2 
dx by 


2 ^15^ * _L 

» ‘ i? 'i i7 


h^ = half thickness of core 


h^ = thickness of facing sheets 


^l’^2 ~ radii of curvature 


= shear modulus of core 


Young's modulus 


density of core 


- density of facing sheets 


Equation (10»1) is based on the following assumptions: 

(1) The facing layers have the same material properties^ are of equal 
thickness, and are much thinner than the core. They carry only 
direct stress and have no flexural rigidity about their own middle 


surfaces. 
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(2) The core layer contributes negligibly to the moment resultants and 
membrane stress resultants of the composite shell. However 
resistance of the core to the transverse shear is quite consider- 
able. 


Now the frequency equation for a sandwich element in terms of the 
wave numbers is given as 


where 


^ ( 


Eh„ 


-) x[ 


, 2 ,, 2 2 . 2 


( Xk^ +k^)^ 

+ ^ ] 




( 10 . 2 ) 


^-R2 


Using the transformation = r cos 0 and = r sin 0 and solving 
(10.2) for Tj it gives 


2 

r 

max 


(oj, 0) 


-Lj (f + [f^ + 4(l-n^)S^f 

2Sh^ 


(10.3) 


where 


fj^( 0 J, 6) = - h^( Xcos^e +sin^e)^/ . 

Here it is assumed that |x| <1 . 

Now the number of resonant modes is given as 

-N(o)) = 5 - J / (u),0)de (10.4) 

2n e^(uj) 
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where the integration is to be carried out for values of 0 over that 
part of the quadrant 0 < 0 < ~ in which the integrand is real and 
positive. 

Hence from equation (10.4) the number of resonant modes is given 
by 


N((ju) 




J (£j + [£j + dS . (10.5) 




Differentiating equation (10.5) with respect to O) gives 


V o2 o2 [f +2(1-^")S"] 

n(m) = — ^ — ( 02 -ei) +— 2—2 ¥ 

2n^Shf UJ ^ i UJ g (uj) [fj +4(l-n^)S^f 


6o(uj) 


.2. „2. 


de 


2 2 *^^2 2 ®1 

+ Sh,[r (u)^ 0„) j - r (u), 0i) j ]] 

I'- max^ ’ 2 dm max^ ^ ^ 


dJT 


(10.6) 


Utilizing the transformation 


2 2 

y = X cos 0 + sin 0 


(10.7) 


equation (10.6) can be written as follows 


2 ^2 


*^1 ^9 o2 

n(uj) = -4-V J 

2TT^Sh? ‘“2 1 (U J 


. .0) ^ 2 c . , 

^ 7^ 

s 2-y 


[(y-X)(l-y)((^)^-y^)(4-y^]^ 


2 2 *^^2 2 

+ Shf[r‘‘ (ou, e,) r (u), 0,) ^1 ) 

I*- max ’ 2 d(u max ^ 1 duj ■' 


“>s 

( 10 . 8 ) 
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where 


and 


2 

m = 7y 

S r, 2 


Eh„ 


[4n (p^lil+p2V^l^ 


2 2^22 
ou = 0) + 0 ) c 

0 s 

2 2 

2 2 ^ 


= X + (1-X)sin^0j^ 

72 = 'X + (l-Xlisin^e^ . 


Equation (10.8) cannot be expressed iramediately in a simpler form 
as it contains hyperelliptical integral. However it can be used to 
obtain the expressions for special cases. 

For a spherical cap 


^1=^2^ X = 1 


and for a flat sandwich plate 


Rl R 2 


Hence the expressions for modal density for a spherical cap and 
flat sandwich plate are given by 


n(u)) = 0 


OD < CD 


s 


'i 'i 2 

n(u)) « (-^)(^U+[p,+2(l-ti^)S^][fJ+4(l-n^)S^fJ"^} 


4rrSh' 


2 JL^2 

for ( 1 ) > tt) . 


(10.9) 
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Figure 10.1 Modal density versus frequency for a sandwich spherical 
cap and flat sandwich plate 
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The variation of modal density versus frequency for spherical cap 
and a flat sandwich plate is illustrated in Figure 4,1 by considering 
the plate and spherical cap having the following material and geometric 
properties: = 60 in. ^ = 36 iUc E = 10 x 10^ PSI, = 0.3, 

h^ = 0.5 in,^ h^ = 0.02 in. = 5.5 Ib/cu ft_, = 170 Ib/cu ft. 


10.3 Orthrotropic Plates 

An orthrotropic plate is characterized by five elastic constants 

E , E , G , LL , Li where E and E are Young's moduli in the x and y 
x^ y^ xy^ xy^ yx x y ° 

direction, G is the shear modulus and u and u represent 
^ xy xy yx 

Poisson's ratio. 

The bending stiffnesses of plate in x and y directions are given 


by 


D 

x 



D 


y 




12hq 


( 10 , 10 ) 


( 10 . 11 ) 


and 

^ n E h^ 

H = (G^yh76) + . (10.12) 

The natural frequencies of free vibrations of a simply supported 
rectangular orthrotropic plate are given by 


(JU 

mn 


. 2 , ^ 
= (in_p!i) 


4 2 2 4 ^ 

[kiD^ + ank^k^ + k^Dy] " 


(10.13) 


where k, and k^ 

1 2 


are the wave numbers. 
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Setting 


1 

= r cos e (10. 14) 

and 

1 

4 

= r sin 9 (10.15) 

then the number of resonant modes are given by 

9 9 (ou) r 

^ ^ 2 ^ ' max 

N(u)) 1/4 X J rdrdG . (10.16) 

TT (Dj^Dy) e^Ctu) 0 

where and are the dimensions of the orthrotropic plate in x and 
y directions respectively. 

The integration of equation (10.6) is carried out over the values 
of 6 in the quadrant 0 < 0 < ^ for which the integrand is real and 
positive. 

Combining equations (10.13)^ (10.14) and (10,15) and solving for 
r gives 


2 

r 

max 


2n(|^l)^ (JD(1- YiSin^2e)"^ 


where 

2y? = 1 - H(D D )"^ . 

^1 X y^ 

For most of the materials 


0 < 


< 


2 ■ 


(10. 17) 


81 



Therefore equation (10.16) reduces to 


1 TT 1 

2D5"5" “2 

N(o))«^(g^) (^) coj (l-YiSin^Ze) d0 


x^c 


0 


(10.18) 


Expressing equation (10.18) in terms of the complete integral of 
the first kind^ it becomes 


i D i 

^ 12 ^ ph ^2 ^ x^4 ^ 

N(cu) (p-^) (g-) F(^ ,, vp 

X c y 


(10.19) 


where 


TT 

2 


(j,k)= f (i-k^sin^e)-^ de . 


Differentiating expression (10.19) with respect to UJ gives 


1,1^ . i D i 

/N 12,ph.2.x. 4 „.tt . 

- — <dV' V *2 ' V 

X c y 


( 10 , 20 ) 


Equation (10.20) gives the expression for modal density of the 
orthrotropic plate. 


10„4 Pretwisted Plates 

A pretwisted plate has a middle surface Z defined by 
Z = ^xy 

where ^ is a pretwist constant and x and y are the Cartesian plate co- 
ordinates. The plate may be considered as hyperbolic paraboloidal 
shell. The plate is shallow if and are small in comparison 

with unity. 
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X 


Figure 10,2 Pretwisted plate 


The natural frequencies of a simply supported shallow pretwisted 
rectangular plate are given by 


(JO 

mn 


2 2 2 2 
Eg % h^Ck^+kj)^ 
= _ c T) r i - 


2 2 2 

2 FT 


( 10 . 21 ) 


where 

k^ and k^ are the wave numbers 
^ is a pretwist constant_, and 
h is a plate thickness. 

Using the transformation 

k^ = r cos 0 and k^ = r sin 0 (10,22) 

it can be found that 
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'L • I (1- 


(10.23) 


where 


2 ^ 

1^1 ,22 

4n pcju 


Hence^ as described in the previous section_, the number of 
resonant modes is given as 


TT 


N((w) 


2VI(l-u 12 ( P J (1- P?sin^2e)^ de P? < 1 . 

TT h ^1 


(10.24) 


Changing equation (10.24) to the standard form of an elliptical 
integral as 


N(u)) 


2V3(l-n^)^V2 p % n 

^ r~ ^ ^^2 ^ h 


(10,25) 


where E(— , k) is a complete elliptical integral of second kind and 

n/ 2 9 2 i- 

is expressed as E(y , k) = J (1 _ k sin 20) d0 . 


Differentiating (10.15) with respect to m, 


n(uj) 


2V3(l-ix^)'^ '^l'^2 p TT 

TT h 4g ^ ^ 


2 * 

U) > (10.26) 

4rr p 


n(u)) =0 


2 

OD^ < (10.27) 

4tt p 
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The expressions (10.26) and (10.27) represent the modal density 
for the pretwisted plate. 

The effect of pretwist on modal density is illustrated in Figure 
10.4 by considering a plate of area 2160 sq in., having the following 
properties: 

£ = 10x10^ PSI, h= 0.25 in., M. = 0.3, p = 170 Ib/ft^. 

10.5 Monocoque Plates under In-Plane Forces 
The natural frequencies of monocoque plates under the action of the 
in-plane forces T^ and (defined as positive in the outward direc- 
tion) are given as follows: 


0) 

mn 


_ 1 ^ 
- 2tt 


2 2 2 2 


12(l-|i ) 





(10.28) 


where 


T 


T 

X 

y 


It is assumed that |T^| ^ |t^| . 

Substituting the transformation = rcos 9 and = rsinG in 
equation (10.18) and solving for r^ it gives 

’^max "" ^ (1 +P2(T cos^ e +sin^0)^]^ 


6(l-U)Tg^ 2 2 

^ — (tcos B+sin'^Q) 


(10.29) 
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Figure 10.4 Modal density versus frequency for a rectangular plate 
under in-plane forces with in-plane forces as parameter 
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where 




2 3(l-/)Tf,g 


y-c 


2 , 2 2 „ 4 

4tt u) pEh 

Again the number of resonant modes can be found as described in 
Section 10.2 and the expression is given as 




TT 

2 


N(0)) 


d0 


c 0 


^ [(1-^ )T (T+1)] . 

4nEh-^ ^ 


(10,30) 


Differentiating (10.30) with respect to oo^ 

H 

n(u)) P. ?yi L i- E ( P )% J [l+p2( Tcos^e+sin^e)^]"^ 


d0. 


c 0 


(10,31) 


Introducing the transformation 


2 2 

y = T cos 0 + sin 0 


the integral can be expressed in terms of the complete elliptic 
integral of the first kind, and consequently the modal density expres- 
sion reduces to the form as shown: 


n(cu) 


2V3(1-h2)% 


TTh 


P 


^Eg '' 


F(J, 


k) 


O 2 +1) 


(10.32) 
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where 


2 _ P2^1-T)^-[(P2+1)^- 

4(P^+1)^ (P^T^+l)^ 

The effects of the in-plane forces on the modal density of a 
rectangular plate is illustrated in Figure 10.3 by considering a plate 
of area 2160 sq in. and having the same material properties as the 
pretwisted plate of the previous section. 


10.6 Monocoque Shells on an Elastic Foundation 
When a rectangular isotropic monocoque shell lies on an elastic 
foundation of modulus its natural frequencies are approximately 
given by 


U) 

mn 


®c 


Eh^(kJ+k^)^ Eh(xk?+k^)^ % 

^ ^ ^ ^ 1 

12(1-/) +kp^ 


(10.33) 


where 



It is assumed that | | < | | or |x| < 1- 

Expression (10.33) is similar to the expression obtained by 
Bolotin for the unsupported shell except the term K. 
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Hence the expression can be written as follows: 


2^ . (1-^^)^ V 2 r p 

TT 7X727¥ h I J 


n(tu) 


62(1“) 


de 


(cu -(JU ) 
s 


0j^(cu) ^3^ 


2 2 ^^2 1 2 2 


n(u)) = 0 


U)^ < cOg (10.35) 


where 


£3(9) = [l-PgCxcos e+sin 6) ]' 


Kg. 


U) 


4TT^ph 


, 2 2,,-l 

E[ — (u) - U)^)] 


Depending on relative magnitudes of X the integral in 

equation (10.34) has different values: 

Considering only positive Gaussian curvature^ there are then 
three subcases, each of which gives different modal densities within a 
certain frequency band. These expressions^ derived by following the 
steps of Bolotin^ are given as: 

* 

For 0) < OD 

s 


n(o)) =0 , 


(10.36) 
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(JU < U) < U) 
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F(j, k) 


Q^/1 V\^ 


(tU^-CDp'* P^(l-X)' 
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* 

U) < (1) 


n(u)) 


2V3(^ V' 


TT 


_1 (_E_)^ X 
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(U 


F( 


r ^ F 


-) 


(0)^.0)^)^ [(l+p3)(l-p3X)]^ 


where 


= O3 +1)(1- P3X)[2P3(1-X)]"^ 


= [X"E( 


. 2.% 


-) + <] 


* 


U) 


4tt^PR? -1 „ % 

[E(— 


The graph of modal density versus frequency is plotted with 

different elastic foundations for a cylindrical panel of area 2160 

6 

sq in. having the following properties: E = 10x10 PSI^ p. = 0.3^ 

h = 0.25 in., p = 170 Ib/cu ft, = 75 in., = <». 


10 . 7 Graphical Results and Discussion 
Figure 10.1 shows that the modal density of a sandwich plate or 
spherical cap has a singularity at the frequency 00 ^. Moreover when the 
shear modulus G of the core is large, the modal density of a plate 
approaches constant, which is the modal density of a monocoque plate 
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r 


2 2 

of bending stiffness [2Ehj^h2 / (l-M' ) ] -^ The modal density of the 
spherical cap approaches asymptotically that of the plate as the 
frequency increases (cju » ou^) . It can also be concluded that as U) 
becomes large^ the modal density of any shallow shell approaches 
as 3 miptotically that of the corresponding plate and above a certain 
frequency^ the modal densities of all sandwich elements increase 
linearly with frequency. 

As seen from expression (10,20) modal density for an orthrotropic 
plate is independent of the frequency (jd^ however it does depend on the 
geometry and material properties of the plate. For a plate having the 
same elastic properties in x and y co-ordinate direction expression for 
modal density is the same as that of the monocoque rectangular plate. 

Figure 10.3 shows the effect of pretwist on the modal density of 
the plate. The pretwist introduces a singularity at the frequency 
* ^ 

(ju = ( ) below which the modal density is zero. Above the 

2tt P 

frequency^ the modal density asymptotically approaches that of an 
untwisted plate_^ which is constant. 

The effect of in-plane forces on a modal density of a plate is 
illustrated in Figure 10,4. It shows that the introduction of in-plane 
forces produces a singularity in the modal density at zero frequency. 
Moreover the modal density of a monocoque plate with in-plane forces is 
greater than the modal density of the unloaded plate^ regardless of the 
sign of forces^ but the amount of increase is not simply related to the 
relative magnitude or sign of the forces. As seen from the graph^ if 
both forces and T^ in x and y co-ordinate directions are equal and 
of opposite sign^ the modal density Is higher than if they were of the 
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same sign. In spite of loading the plate^ as the frequency increases^ 
the modal density of the loaded plate approaches asymptotically that 
of the unloaded plate. 

For the monocoque shells of positive Gaussian curvature^ placed 

on an elastic foundation^ variation of the modal density is illustrated 

in Figure 10.5. It shows that the modal density is zero below the 
* 

frequency cu at which it has a singularity. The modal density has 
s 

•/r j- 

second singularity at the frequency cw . However above the modal 

density of the shell decreases monotonically and asymptotically 

approaches that of the corresponding plate as oj becomes very large. 

The foundation modulus k modifies the position of the singularities 

according to relative magnitudes of uj* and cju . 

s 
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11c SUMMARY AND CONCLUSIONS 


The modal density of a structural element is defined as the number 
of resonant modes within a unit frequency intervale In this thesis^ 
the expressions and graphs are presented which can be used to estimate 
the average modal densities of structural elements like rods^ beams^ 
plates^ thin cylindrical^ spherical and conical shells^ composite 
structures^ shallow sandwich shells^ orthrotropic plates^ pretwisted 
plates, plates subject to in~plane forces and shells resting on elastic 
foundation^ The expressions and graphs are valid for elements of 
arbitrary shape and having any prescribed boundary conditions. 

In case of circular rods having uniform cross sections, modal 
density is independent of geometry of rod and frequency of vibrations 
for both longitudinal as well as torsional vibrations. However, modal 
density for torsional vibrations is about 1„ 5 times that of longi- 
tudinal vibrations . 

The modal density of beams undergoing transverse vibrations 
depends both on geometry of beam and frequency of vibration. It 
decreases with the frequency and reaching asymptotically to zero value 
for large frequencies. 

In case of solid flat rectangular and circular plates the modal 
density is constant for a given plate. Thus it is independent of 
frequency of vibration but not the geometry of the plate. For a given 
frequency, modal density of rectangular plate is approximately equal to 
that of the circular plate, both the plates having equal area, and 
thickness and of the same material. 
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Modal density of thin cylindrical shells decreases above the ring 
frequency^ reaching a constant value asymptoticallye However^ below 
ring frequency there is a variation that is linear on a log-log ploto 
The modal density for thin spherical shells has a singularity at 
the ring frequency^ below which the modal density is zero. However, 
above ring frequency it decreases monotonically, approaching that of a 
plate asymptotically as frequency becomes large. 

The thin conical shell shows that expressions presented for modal 
density are applicable over a wide range of the cone geometries and 
frequency ranges of practical interest. 

The modal density of a composite structure^ as derived analytical- 
ly^ is additive over its components. 

For the shallow sandwich elements like the spherical cap modal 
density increases with increasing frequency^ and for very large values 
of frequency modal density of the sandwich shells approaches asymp- 
totically that of sandwich plates. 

In case of the orthrotropic plates^ modal density is frequency 
independent but it does depend on the geometry and material properties. 
The modal density of pretwisted plate decreases monotonically 
after some value of frequency and asymptotically approaches that of 
the untwisted plate as frequency becomes large. By increasing the 
pretwist constant the value of frequency at which singularity occurs 
also increases. 

The modal density of a monocoque plate with in-plane forces is 
greater than the modal density of an unloaded plate^ regardless of the 
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sign of forces and as the frequency increases^ the modal density of 
the loaded plate approaches asymptotically that of the unloaded plate„ 

The modal density of shells on an elastic foundation decreases 
monotonically and asymptotically approaches that of the corresponding 
plate as frequency becomes large„ 

Thus in general modal density of all elements except plates and 
rods are frequency dependent. The modal density of the shells 
approaches asymptotically that of the corresponding plate. For the 
sandwich elements the modal density increases with increasing 
frequency^ whereas the modal density of monocoque elements approaches 
a constant value. 

The concept of modal density is very useful in solving the multi- 
modal vibration problems and is of great value when the input force 
or excitation is random. In this type of analysis without the apparent 
knowledge of mode shapes and frequencies it is possible to give some 
insight into the response of the structures to the given excitation and 
some insight into the amount of energy which will be absorbed. More- 
over modal density of structures is relatively independent of the 
boundary condition; it is a useful tool in estimating average response 
levels of multimodal vibration. 
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13. APPENDIX. LIST OF SYMBOLS 


a 




a 

n 


D ,D 

y 


E .E 

y 

G 

h 

h. 


I 

J 

r 

K 

k 


k k 




radius of cylinder or circular plate^ in. 
principal dimensions of shell surface^ in. 
nn/ 


longitudinal wave velocity = 

torsional wave velocity = 

3 2 

stiffness modulus = Eh / (12(l-|i. ) 

stiffness moduli in x and y co-ordinate direction 

Young*s modulus of elasticity 

Young's moduli in x and y co-ordinate direction 

shear modulus of elasticity 

gravitational constant 

thickness of shell wall or plate 

half thickness of core in sandwich element 

thickness of facing sheets in sandwich element 

moment of inertia of a section 

Bessel function of first kind of order n 

modulus of elastic foundation 

constant 

wave number 

length of cylinder or beam 

surface dimensions of the plate 

length of cone^ apex to base slant length 
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m 


N(cu) 

N(V) 

n(ou) 

n(v) 

n 


^0 

R2 

r 


T ,T 

y 

V 

V 


Xi,X2 

X(x) 

Y 

n 


length of cone truncation^ apex to top slant length 
integer value or number of circumferential waves 
number of resonant modes 

number of resonant modes in terms of dimensionless frequency 
modal density 

modal density in terms of dimensionless frequency 

integer value or number of one-half longitudinal waves 

one-half circumferential modes 

amplitude of excitation 

radius of shell curvature 

cylindrical co-ordinate 



time 

in-plane forces in x and y co-ordinate direction 

radial velocity amplitude 

axial velocity amplitude 

tangential velocity amplitude 

displacement normal to surface 

generalized co-ordinates 

normal mode function 

rectangular co-ordinates 

Bessel function of second kind of order n 


stress coefficient 
truncation ratio = L^/L 
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p 


Ak 

Ak 

e 


1 

2 


K 


V 


h/2V3a 

change in longitudinal wave number 
change in circumferential wave number 
cylindrical co-ordinate 
radius of gyration 
Poisson’s ratio 

dimensionless frequency - ^ 


3.14 


P 

PpPz 

O’ 

T 

X 

t 

U) 


density of material 

density of core and facing sheets 

mrra/ ^ 

T 

X 

t“ 

y 

stress function or angular displacement 

Rj/R^ 

one-half cone angle at apex 
angular frequency 


F(* 2 - ; k) complete elliptic integral of first kind 

E(~ j k) complete elliptic integral of second kind 


NASA-Langley, 1971 32 CR-1775 
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